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HIGHLIGHTS 


►  The  developed  mathematical  model  can  predict  direct  methanol  fuel  cell  performance. 

►  The  model  is  validated  using  data  at  various  operating  conditions. 

►  This  model  can  be  extended  to  methanol  fuel  cells  with  various  membranes. 

►  The  model  supplies  a  cost-  and  time-saving  advantages  for  researchers. 
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One-dimensional  mathematical  model 


The  objective  of  this  study  is  to  develop  a  one-dimensional  mathematical  model  for  direct  methanol  fuel 
cell  (DMFC)  performance  prediction.  The  modified  model  can  be  applied  to  DMFCs  operated  at  wider 
ranges  of  methanol  concentrations  and  temperatures  with  various  solid  electrolytes.  Based  on  the  model 
of  Li  et  al.  [J.  Power  Sources  154  (2006)  115-123],  we  propose  a  procedure  to  predict  the  effective 
diffusion  coefficients  of  reactants  (methanol  and  oxygen)  in  the  diffusion  layers  to  further  estimate  the 
limiting  current  densities.  This  approach  does  not  use  empirical  fitting  parameters  and  adopts  the 
Darken  equation  and  Boltzmann  relationship  to  express  the  diffusivity  dependency  on  the  anode  and 
cathode  feed  compositions.  The  diffusion  layer  porosity  effect  on  the  methanol  and  oxygen  diffusion 
coefficients  is  accounted  for  using  the  Das  equation.  The  proposed  model  is  validated  by  experimental 
results  obtained  for  1-5  M  of  methanol  feed  at  operating  temperatures  of  30-90  °C  using  Nation, 
sulfonated  poly(phthalazinone  ether  ketone),  and  sulfonated  poly(ether  ether  ketone)  electrolytes.  The 
predictions  are  in  good  agreement  with  the  experimental  data.  This  model  is  helpful  for  researchers  in 
estimating  cell  performance  by  eliminating  the  costs  involved  with  preparing  single  cells  during  new 
material  development. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Direct  methanol  fuel  cells  (DMFCs)  are  a  promising  alternative 
power  source,  especially  for  portable  electronic  equipment,  such  as 
laptop  computers,  two-way  radios,  and  mobile  phones.  DMFCs 
have  the  advantages  of  high  energy  density,  low  emission  of 
pollutants,  and  the  ability  to  use  liquid  fuel,  which  is  easier  to  store 
and  transport  than  hydrogen  [1-4].  Moreover,  fuel  cells  can  be  used 
as  an  emergency  power  supply  because  of  their  rapid  start-up  at 
ambient  temperatures. 


*  Corresponding  author.  Tel.:  +886  3  2118800x5489;  fax:  +886  3  2118700. 

E-mail  address:  jessie@mail.cgu.edu.tw  (S.J.  Lue). 

1  Equal  contribution  to  first  author. 

0378-7753 /$  —  see  front  matter  ©  2012  Elsevier  B.V.  All  rights  reserved. 
http://dx.doi.org/10.1016/jjpowsour.2012.ll.029 


In  the  past  decades,  substantial  research  has  been  devoted  to 
establishing  models  to  better  understand  the  phenomena  arising 
within  DMFCs.  Most  models  consider  the  mass  transport  and 
electrochemical  properties  in  a  DMFC.  Sundmacher  and  Scott  re¬ 
ported  on  a  one-dimensional  model  for  DMFCs  using  a  vaporized 
methanol  feed  [5].  Scott  et  al.  [6-8]  also  developed  a  one¬ 
dimensional  model  for  liquid-feed  DMFCs  that  accounts  for  the 
methanol  crossover  and  mass-transfer  limitation  at  the  catalyst 
surface.  The  influence  of  methanol  crossover  from  the  anode  to 
cathode  is  based  on  a  combination  of  diffusion,  electroosmotic 
drag,  and  pressure.  Kulikovsky  [9,10]  described  a  one-dimensional 
semi-empirical  model  for  DMFC  cell  performance  based  on  Tafel 
kinetics  of  the  electrochemical  reaction  of  methanol  oxidation,  the 
diffusion  transport  of  methanol  through  the  backing  layer,  and  the 
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methanol  crossover.  These  models  have  several  (5-9)  fitting 
parameters  and  were  solved  under  a  narrow  range  of  methanol 
feed  concentrations  (0-0.75  M).  Eccarius  et  al.  [11  ]  used  a  design  of 
experiments  coupled  with  a  mathematical  model  for  a  DMFC  to 
study  the  influence  of  temperature,  membrane  thickness,  methanol 
molarity,  anode  flow  rate  on  methanol  crossover,  and  leakage 
current  density.  Ko  et  al.  [12]  studied  the  numerical  simulation 
results  of  a  DMFC  model  and  compared  them  with  experimental 
results  on  the  methanol  crossover  density  of  different  Nation 
membranes.  Jeng  and  Chen  [13]  further  defined  the  methanol 
concentration  profile  by  introducing  the  void  fractions  of  the 
diffusion  and  catalyst  layers  and  the  void  fractions  of  Pt-Ru/C  and 
ionomers  in  the  catalyst  layer.  Kareemulla  and  Jayanti  [14]  used  the 
multi-step  reaction  mechanism  to  describe  the  electrochemical 
oxidation  of  methanol  at  the  anode  and  the  Stefan-Maxwell 
equations  for  multi-component  diffusion  on  the  cathode  side. 
Sun  et  al.  [15]  and  Ge  and  Liu  [16]  introduced  a  three-dimensional, 
two-phase  transport  model  for  liquid-feed  DMFCs.  A  review  of  the 
literature  on  DMFC  models  is  presented  by  Oliveira  et  al.  [17]. 

All  the  above  models  predict  DMFC  cell  performance  at  low 
methanol  concentrations,  and  most  of  the  studies  compare  the 
modeling  results  with  fuel  cells  using  Nation  electrolyte.  Some  of 
the  transport  parameters  are  limited  in  the  literature. 

In  this  paper,  we  focused  on  developing  a  mathematical  model 
to  expand  the  applicability  of  DMFC  performance  prediction.  The 
proposed  model  is  based  on  the  Kulikovsky  model  [9],  which  was 
later  detailed  by  Li  et  al.  [18];  however,  our  model  incorporates 
analytical  mass  transfer  parameters  to  precisely  estimate  the 
effective  diffusion  coefficients  of  methanol  and  oxygen  in  the 
diffusion  layers.  A  benefit  of  this  approach  is  that  the  reactant 
diffusion  coefficient  at  certain  feed  compositions  and  cell  temper¬ 
atures  can  be  used  to  predict  the  diffusivity  under  other  operating 
conditions.  This  model  takes  into  consideration  various  methanol 
concentrations,  temperatures,  and  different  solid  electrolytes. 
Furthermore,  the  DMFC  performance  obtained  by  our  model  is 
validated  using  experimental  data. 

2.  Model  development 


(i)  The  cell  is  under  isothermal  conditions. 

(ii)  The  methanol  concentration  in  the  anode  catalyst  layer  is 
constant. 

(iii)  The  electrode  kinetics  can  be  described  using  the  Tafel 
equation. 

(iv)  Methanol  crossover  is  mainly  caused  by  diffusion  in  the 
membrane  electrolyte;  the  electro-osmosis  effect  on  the 
methanol  transport  rate  through  the  membrane  is  negligible. 

(v)  The  methanol  concentration  gradient  is  linear  with  respect  to 
the  trans-membrane  direction  inside  the  membrane. 

(vi)  The  gas  hold-up  volume  due  to  carbon  dioxide  evolution  is 
negligible;  the  produced  gas  is  immediately  removed  by  the 
anode  feed  flow. 

(vii)  Membrane  electrolyte  resistance  governs  the  ohmic  voltage 
drop  in  the  cell. 

(viii)The  oxygen  pressure  in  the  cathode  catalyst  layer  is  constant. 

(ix)  The  reactions  at  the  anode  and  cathode  are  first  order. 

(x)  The  concentration  distributions  are  one-dimensional  and  only 
vary  with  respect  to  the  trans-membrane  direction. 

The  overpotential  of  a  DMFC  derived  from  Kulikovsky  [9]  and  Li 
et  al.  [18]  is  expressed  by  the  following  equation: 

v  =  E  -  77a  -  Vti  -  Vc  (2) 

where  F  is  the  thermodynamic  potential  of  a  DMFC  at  a  given 
operating  condition,  is  the  overpotential  at  the  anode  described 
in  Section  2.1.1,  ric  is  the  overpotential  at  the  cathode  described  in 
Section  2.1.2,  and  t/q  is  the  ohmic  potential  drop  in  the  cell.  7]q  is 
calculated  as  [9]: 

iL  f  . 

Vn  =  -  (3) 

where  i  is  the  current  density,  L  is  the  membrane  electrolyte 
thickness,  and  a  is  the  ionic  conductivity  of  the  electrolyte. 
According  to  the  Nernst  equation,  the  thermodynamic  cell  potential 
at  any  given  condition  (F)  is: 


Our  proposed  mathematical  model  is  based  on  the  models 
proposed  by  Kulikovsky  [9]  and  Li  et  al.  [18]  but  uses  mass-transfer 
characteristics  from  analytical  equations.  First,  we  adopted  the 
Wilke-Chang  estimation  method  [19],  the  Lennard-Jones  potential 
model  [20],  and  the  equation  of  state  of  vapor-liquid  equilibrium 
[21]  to  estimate  the  methanol  and  oxygen  diffusion  coefficients  in 
the  fluids.  These  diffusion  coefficients  were  modified  in  accordance 
with  the  formula  obtained  by  Das  et  al.  [22]  to  calculate  the 
effective  diffusion  coefficients  of  methanol  and  oxygen  in  the 
porous  diffusion  electrodes.  The  diffusion  coefficients  were  used  to 
predict  the  limiting  current  densities  under  various  operating 
conditions.  Finally,  we  applied  the  Kulikovsky  [9]  model  to  predict 
the  cell  voltages  and  to  calculate  the  power  densities  in  a  DMFC 
using  a  1  -5  M  methanol  feed,  operating  temperatures  ranging  from 
30  °C  to  90  °C,  and  Nation  and  other  sulfonated  solid  electrolytes. 
The  developed  model  is  described  in  the  following  sections. 

2.2.  Conventional  model 


E 


E° 


RT  .  QC02aH20 
nF  3 

aCH3OHao2 


(4) 


where  a  is  the  activity  of  each  reactant  or  product,  R  is  the  gas 
constant,  T  is  the  absolute  temperature  in  Kelvin,  n  is  the  number  of 
electrons  in  the  half-reaction  in  a  DMFC  (6  in  Eq.  (1)),  Fis  Faraday’s 
constant,  and  F°  is  the  thermodynamic  cell  potential  of  the  DMFC 
under  standard  conditions  (a  pure  substance  at  25  °C  and  1  bar): 


F°  = 


A  G° 
nF 


(5) 


where  AG°  is  the  Gibbs  free  energy  change  of  the  methanol 
oxidation  reaction  in  Eq.  (1).  The  predicted  power  density  (P)  of 
a  DMFC  was  calculated  as  the  product  of  the  predetermined  current 
density  and  the  estimated  cell  voltage  (from  Eq.  (2)): 


The  overall  reaction  of  a  DMFC  is: 


P  =  iV 


(6) 


CH3OH(1)  +  |o2(g)^C02(g)  +2H20(1)  (1) 

The  assumptions  for  this  model  are  described  in  the  following 
sections: 


2.2.2.  Anode  overpotential 

The  anode  overpotential  (r/a)  can  be  calculated  by  the  combi¬ 
nation  of  the  kinetic  reaction,  mass  transport,  and  methanol 
crossover  as  described  in  Eq.  (7)  [18,23]: 
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Va 


jba  In 


bi  lnl  1 


+  b a  ln(l  +  fi) 


(7) 


(15) 


where  ba  is  the  Tafel  slope  at  the  anode,  i  is  the  experimental 
current  density,  iea  is  the  exchange  current  density  at  the  anode, 
and  i'ia  is  the  limiting  current  density  of  the  anode  and  is  defined 
as  [24]: 

I'la  =  6F^ba^ni  (8) 

^ba 

where  Dba  is  the  effective  diffusion  coefficient  of  methanol  aqueous 
solution,  Lba  is  the  thickness  of  anode  diffusion  layer,  and  Cm  is  the 
methanol  feed  concentration. 

In  Eq.  (7),  ^  is  a  dimensionless  parameter  and  is  calculated  as: 


and 


iec  =  0.04222 


/73, 200' 
eXP  014 


1 


273  +  80 


(16) 


This  relationship  is  based  on  an  activation  energy  of 
73,200  J  mol-1  for  a  DMFC,  conducted  between  30  and  80  °C  [26] 
using  Nafion-117  sandwiched  between  Pt  microelectrodes.  The 
membrane  characteristics  play  important  roles  in  the  cell  perfor¬ 
mance.  For  instance,  electrolyte  conductivity  (er),  methanol 
permeability  (0),  and  membrane  thickness  (L)  are  shown  in  Eqs.  (3) 
and  (9).  The  data  are  further  entered  to  Eqs.  (2),  (12)  and  (14)  to 
demonstrate  their  impacts  on  the  fuel  cell  voltage. 


lb* 

L  ^ba 


(9) 


where  (3  is  the  methanol  permeability  of  the  membrane  electrolyte. 
Wang  and  Wang  demonstrated  that  the  ba  and  iea  values  are 
functions  of  temperature  [24]: 


2.2.  Modified  model 

The  equations  stated  in  Section  2.1  were  adopted  in  the  modi¬ 
fied  model  for  DMFC  performance  prediction  with  parameter 
adjustments.  The  solid  electrolyte  conductivity  (a)  as  a  function  of 
temperature  (T)  is  described  as  follows: 


ba 


RT 
08 F 


and 


lea 


94.25 


/35, 570 
exp( 


1 


8.314  J  1273  +  80 


(17) 

where  crref  is  the  conductivity  at  a  reference  temperature  (Tref)  [27]. 
The  Dba  and  Dbc  values  were  adjusted  using  the  procedure 
described  in  the  following  sections. 


(10) 


a  =  (T ref  exp 


1268 


1  1 


Tref  T 


where  35,570  (in  J  mol-1)  is  the  activation  energy  derived  by  linear 
regression  from  the  data  of  the  performance  of  a  DMFC  that  was 
composed  of  a  3.0  mg  cm-2  Pt-Ru  anode,  a  3.4  mg  cnrT2  Pt  cathode, 
and  a  Nafion-117  electrolyte.  The  value  of  8.314  (in  J  mol1  K-1) 
represents  the  gas  constant.  These  relationships  were  established 
at  60-100  °C  using  1  M  methanol  (4  mL  min-1)  and  2.07  x  105  Pa 
(30  psig)  air  at  1800  seem  [25]. 

2.2.2.  Cathode  overpotential 

The  cathode  overpotential  (r]c)  caused  by  the  diffusional  mass- 
transfer  limitation  due  to  a  loss  in  oxygen  concentration  can  be 
estimated  from  the  combination  of  the  kinetic  reaction,  mass 
transport,  and  methanol  crossover  using  the  following  equation: 

r?c  =  bcln(4-)  -bclnfl  -J--Rc)  (12) 

where  bc  is  the  Tafel  slope  at  the  cathode,  iec  is  the  exchange  current 
density  of  the  cathode,  i\c  is  the  limiting  current  density  at  the 
cathode,  and  Rc  is  a  dimensionless  parameter  accounting  for  the 
potential  drop  due  to  the  methanol  crossover  [24].  The  tjc  and  Rc 
values  were  obtained  using  the  following  equations: 

he  =  (13) 

^bc 


2.2.2.  Methanol  diffusion  coefficient  at  the  anode 

The  methanol  diffusion  coefficient  in  aqueous  solution  was 
estimated  using  the  Darken  equation  [28,29]: 

^MW  —  (^MW^M  +^MW%)ai  (IS) 

where  Dmw  is  the  methanol  diffusion  coefficient  in  water  from 
methanol  solutions,  DfJ  is  the  diffusion  coefficient  of  solute  /  in  an 
infinite  dilute  solution  of  solvent  J ,  and  M  and  W  denote  methanol 
and  water,  respectively.  The  xm  and  Xw  values  are  the  methanol  and 
water  molar  fractions,  respectively.  The  parameter  a  is  a  thermo¬ 
dynamic  correction  term  given  by  the  following  equation  [30]: 

« =  (i9) 

\d  In  Xm/  j  p 

where  om  and  xm  are  the  activity  and  mole  fraction  of  methanol, 
respectively.  The  D^M  and  D^w  values  were  predicted  using  the 
Wilke-Chang  equation  for  the  diffusion  coefficients  of  solutes  in 
solvents  at  infinite  dilution  [19]: 


n° 


7.4  x  10 -8(<2>mMm)1/2T 

?m46 


and 


(20) 


where  Dbc  is  the  diffusion  coefficient  of  oxygen  in  the  cathode 
backing  layer,  Lbc  is  the  thickness  of  the  cathode  backing  layer  and 
Co  is  the  oxygen  feed  concentration,  and 


n° 

^MW 


7.4  x  10-8(<FwMw)1/2r 
Cwi'm6 


(21) 


Rc 


ila 

he 


(14) 


The  bc  and  iec  values  depend  on  the  cell  temperature  [24]: 


where  M  is  the  solvent  molecular  weight,  v  is  the  solute  molar 
volume  at  its  boiling  temperature,  J  is  the  solvent  viscosity,  T  is  the 
temperature,  and  $  is  the  association  factor  of  the  solvent  (2.6  for 
water  and  1.9  for  methanol  [19]). 
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The  correction  term  a  in  Eq.  (19)  was  acquired  from  the  vapor- 
liquid  equilibrium  data  of  methanol  and  water  [21  ]  using  UNIQUAC 
parameters.  The  In  cjm  and  In  xm  from  30  to  80  °C  are  plotted  in 
Fig.  1.  From  the  slope,  we  obtained  a  values  of  0.9704,  0.9419, 
0.9146,  and  0.8639  for  1  M,  2  M,  3  M,  and  5  M  methanol  solutions, 
respectively,  at  30  °C.  As  the  temperature  increased,  the  a  values 
slightly  decreased.  For  instance,  the  a  values  were  0.9594,  0.9211, 
0.8851,  and  0.8206  for  1  M,  2  M,  3  M,  and  5  M  methanol  solutions, 
respectively,  at  80  °C. 

By  substituting  appropriate  methanol  concentrations  and 
temperatures,  we  obtained  corresponding  methanol  diffusion 
coefficients  in  water.  Then,  the  effective  methanol  diffusion  coef¬ 
ficient  in  the  porous  diffusion  layer  was  calculated  as  [22,31]: 

<w  =  (l  -  CttOK  (22) 


(25) 


where  Mo  and  Ma  are  the  molar  masses  of  oxygen  and  air, 
respectively.  is  the  diffusion  collision  integral  [32]: 


1.06036  0.193  1.03587 

(7*)°i56i  +exp(0.47635r)  +exp(l. 529967*) 

1.76474 

+  exp(3.8941lT) 


(26) 


where  T  (reduced  temperature)  equals  kT/e oa,  k  is  Boltzmann’s 
constant,  and  eoa  is  the  characteristic  Lennard-Jones  energy  pre¬ 
sented  in  the  following  equation: 

£oa  =  (^A)1/2  (27) 


where  6a  is  the  porosity  of  the  anode  diffusion  layer.  The  experi¬ 
mental  methanol  diffusion  coefficient  at  90  °C  has  been  reported  to 
be  1.8  x  10-5  cm2  s_1  with  a  ba  of  0.3  mm  [9].  This  value  was  used  to 
calculate  the  apparent  porosity  0a  (which  was  estimated  to  be 
0.406).  The  predicted  effective  methanol  diffusion  coefficient  was 
used  as  Dba  in  Eqs.  (8)  and  (9)  to  estimate  the  anode  overpotential 
described  in  Section  2.1.1. 


where  ro  and  ea  are  the  characteristic  Lennard-Jones  energies  of 
oxygen  and  air  [30]. 

The  effective  oxygen  diffusion  coefficient  (D^)  in  the  porous  gas 
diffusion  layer  was  calculated  as  [22]: 

<  -  (’  -  (t^))d“  <28> 


2.2.2.  Oxygen  diffusion  coefficient  into  cathode 

The  oxygen  diffusion  coefficient  in  air  at  the  cathode  was  esti¬ 
mated  using  the  simplified  Boltzmann  equation  [30]: 


0.00266T15 


(23) 


where  T  is  the  temperature,  p  is  the  oxygen  pressure,  Moa  is  the 
equivalent  molar  mass  (defined  in  Eq.  (25)),  Qo  is  the  dimensionless 
diffusion  collision  integral  presented  in  Eq.  (26)),  and  <7oa  is  the 

o 

characteristic  Lennard-Jones  length  (A)  and  is  taken  as  the  average 
of  oxygen  and  air  data: 


where  6C  is  the  porosity  of  the  cathode  gas  diffusion  layer.  The 
experimental  effective  diffusion  coefficient  of  oxidant  at  the 
cathode  at  90  °C  has  been  reported  to  be  0.9  x  1CT3  cm2  s-1  [10] 
when  Lbc  is  0.3  mm.  This  value  was  used  to  calculate  the 
apparent  porosity  0C.  Once  the  apparent  porosity  is  known,  the 
effective  gas  diffusion  coefficient  in  the  porous  electrode  can  be 
determined  and  used  to  estimate  the  cathode  overpotential 
(described  in  Section  2.1.2). 

3.  Experimental 

3.1.  Electrolyte  preparation 


aOA 


<Jq  +  crA 


(24) 


where  gq  and  g a  data  were  obtained  from  the  literature  [28].  Moa  is 
the  equivalent  molar  mass,  defined  as  [30]: 


In  X 


The  Nafion-117  membrane  (from  DuPont  Co.,  Fayetteville,  North 
Carolina,  USA)  was  first  boiled  with  hydrogen  peroxide  (5  wt%, 
Riedel-de  Haen,  Seelze,  Germany)  for  1  h  to  remove  organic 
matters.  It  was  rinsed  in  pure  water  for  0.5  h.  Then,  the  Nafion  was 
boiled  in  H2S04  (1  M,  Riedel-de  Haen,  Seelze,  Germany)  for  1  h  to 
ensure  complete  conversion  to  an  H-type  exchange  membrane  and 
rinsed  again  with  pure  water.  The  pure  water  was  produced  using 
a  Millipore  water  purifier  (RiOs-5  and  Millipore-Q.  Gradient,  Milli- 
pore  Corp.,  Bedford,  Massachusetts,  USA). 

The  sulfonated  poly(phthalazinone  ether  ketone)  (SPPEK) 
powder  with  an  ion-exchange  capacity  of  1.42  meq  g'1  was  ob¬ 
tained  from  FuMA-Tech  GmbH  (Fumion®  P-700,  Ingbert,  Germany). 
In  total,  1  g  of  SPPEK  was  dissolved  in  9  g  of  dimethylformamide 
(DMF)  with  continuous  stirring  for  one  day.  The  solution  was 
filtered  through  Advantec  filter  paper  (No.  1,  Advantec  MFS,  Inc., 
Dublin,  CA,  USA)  to  obtain  a  clear  solution.  The  polymer  solution 
was  cast  on  a  glass  plate  and  the  solvent  was  allowed  to  evaporate 
in  a  ventilation  hood  for  3  days.  The  film  was  then  dried  at  80  °C  in 
a  vacuum  oven  for  24  h.  The  dry  film  was  peeled  off  and  immersed 
in  1  M  H2SO4  (Riedel-de  Haen,  Seelze,  Germany)  for  1  h  to  ensure 
complete  conversion  to  H-type  exchange  membrane  and  was 
rinsed  again  with  pure  water. 

3.2.  Conductivity  and  methanol  permeability  measurements 


Fig.  1.  Plots  of  in  aM  compared  with  in  xM  for  methanol  solutions  at  30-80  °c.  The  The  tested  membrane  electrolyte  was  immersed  in  deionized 

slopes  at  any  methanol  mole  fraction  correspond  to  the  a  value  defined  in  Eq.  (15).  water  or  methanol  solutions  for  3  h  at  various  temperatures  and 
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held  between  two  stainless  steel  electrode  rods  in  a  glass  T-tube,  as 
described  in  our  previous  paper  [33].  The  specimen  was  maintained 
at  95%  relative  humidity  at  a  predetermined  temperature  for  the 
conductivity  measurement.  An  alternate-current  impedance  spec¬ 
trum  was  obtained  using  a  potentiostat  (Autolab  PGSTAT  30ZN 
potentiostat,  Eco  Chemie  B.V.,  Utrecht,  The  Netherlands).  The 
impedance  was  scanned  at  106  Hz-ICT1  Hz  at  an  applied  voltage  of 
10  mV.  Nyquist  plots  were  used  to  extract  the  conductivity  values  of 
the  electrolytes,  as  described  in  our  previous  report  [33]. 

The  methanol  permeability  was  determined  using  a  hand-made, 
double-jacked,  glass  permeation  cell.  The  glass  cell  was  divided  into 
two  compartments.  One  reservoir  was  filled  with  a  1  M  methanol 
aqueous  solution,  and  the  other  receiving  reservoir  was  filled  with 
deionized  water.  The  membrane  was  placed  between  these  two 
compartments  and  the  methanol  transport  through  the  membrane 
was  measured  by  the  concentration  increase  in  the  receiving 
reservoir  with  time.  The  methanol  concentration  was  measured 
using  a  gas  chromatograph  (HP  4890A,  Agilent  Technologies  Co. 
Ltd.,  St.  Louis,  MO,  USA).  The  experimental  set-up  and  calculation  of 
the  methanol  permeability  were  described  in  our  previous  publi¬ 
cation  [34]. 

3.3.  Fuel  cell  test 

The  membrane  electrolyte  was  sandwiched  between  the  sheets 
of  the  gas  diffusion  electrodes  (E-tek,  6  mg  cm-2  Pt-Ru  for  the 
anode  and  5  mg  cm-2  Pt  for  the  cathode)  to  obtain  a  membrane 
electrode  assembly  (MEA).  The  active  area  of  the  MEA  electrode 
was  5  cm2.  Two  heating  tapes  were  adhered  to  the  surfaces  of  the 
end  plates.  The  thermostated  methanol  solution  was  fed  into  the 
anode  at  a  flow  rate  of  2  mL  min-1,  and  humidified  oxygen  gas  was 
fed  into  the  cathode  at  a  flow  rate  of  80  mL  min-1.  A  thermocouple 
was  inserted  into  the  end  plates,  and  the  measured  temperature 
was  fed  back  to  a  controller  to  ensure  that  the  cell  temperature  was 
at  the  set  point.  The  electrochemical  measurements  of  the  DMFCs 
were  recorded  using  the  constant  current  density  mode  with  an 
electrical  load  (PLZ164  WA  electrochemical  system,  Kikusui  Elec¬ 
tronics  Corporation,  Tokyo,  Japan).  The  experimental  set-up  of  the 
fuel  cell  test  was  reported  previously  [35]. 

4.  Results  and  discussion 

4.1.  Effect  of  diffusion  layer  porosity  on  cell  performance 

The  diffusion  layer  porosity  has  a  profound  impact  on  the 
effective  reactant  diffusion  coefficients  in  the  diffusion  layers  (Eqs. 
(22)  and  (28))  and  on  the  cell  voltage  (Eq.  (2)).  Fig.  2  plots  the 
predicted  cell  voltage  as  a  function  of  the  current  density  with  the 
same  porosity  value  for  the  anode  and  cathode  diffusion  layers 
using  the  modified  DMFC  model,  which  was  fed  with  1  M  methanol 
and  1-bar  pure  oxygen  at  60  °C.  When  the  porosities  were  both  0.1 
for  the  anode  and  cathode  diffusion  layers,  the  predicted  voltage 
values  were  the  lowest  among  the  tested  parameters. 

As  the  diffusion  layer  porosities  increased  beyond  0.4,  the  cell 
voltages  increased,  in  particular  for  current  densities  beyond 
70  mA  cm-2  (Fig.  2).  Moreover,  the  polarization  curves  with 
porosities  of  0.4-0.9  were  almost  undifferentiated,  which  implies 
that  the  predicted  cell  voltage  was  nearly  independent  of  porosity 
for  values  higher  than  0.4.  The  predicted  polarization  curves  were 
in  accordance  with  the  experimental  data  (Fig.  2),  except  for  the 
low  current  density  (<50  mA  cm-2)  range.  In  the  literature 
[6,7,12,13,16,24,36,37],  porosity  values  of  0.3-0.875  were  adopted 
in  establishing  models  for  predicting  DMFC  performance;  the 
selection  of  this  parameter  value  only  plays  a  minor  role  based  on 
our  results  presented  in  Fig.  2. 


Fig.  2.  Predicted  voltage-current  density  curves  using  the  modified  model  with  the 
same  porosities  (0.1 -0.9)  for  anode  and  cathode  in  a  DMFC  fed  with  1  M  methanol  and 
pure  oxygen  at  60  °C.  The  experimental  data  (Exp.)  are  shown  for  comparison. 

Next,  we  set  the  anode  and  cathode  porosities  to  different 
values  to  examine  their  individual  impacts  on  cell  performance. 
The  modeling  was  performed  for  a  DMFC  fed  with  1  M  methanol 
and  pure  1-bar  oxygen  at  60  °C.  When  the  cathode  porosity  was 
set  at  0.12  and  the  anode  porosity  varied  from  0.1  to  0.9,  the 
predicted  cell  voltage  was  higher  with  a  higher  anode  porosity 
(Fig.  3(a)).  The  estimated  peak  power  density  (Pmax)  increased 
from  41.7  to  74.5  mW  cm-2  as  the  anode  porosity  increased  from 
0.1  to  0.9.  When  the  cathode  porosity  was  set  at  0.41  and  the 
anode  porosity  varied  from  0.1  to  0.9,  the  predicted  cell  voltage 
also  increased  with  a  higher  anode  porosity  (Fig.  3(b)).  These 
polarization  curves  were  similar  to  those  in  Fig.  3(a)  at  the  same 
anode  porosity  value,  and  the  estimated  Pmax  was  less  than  0.5% 
deviation  from  the  results  using  a  cathode  porosity  of  0.12.  This 
finding  indicates  that  the  cathode  porosity  is  not  a  governing 
factor  in  cell  performance. 

With  the  anode  porosity  set  at  0.41  and  0.12,  the  predicted 
polarization  curves  are  illustrated  in  Fig.  4(a)  and  (b),  respectively. 
At  an  anode  porosity  of  0.41,  the  generated  current  density  reached 
at  least  300  mA  cm-2  (Fig.  4(a))  with  an  estimated  Pmax  of 
70  mW  cm-2,  regardless  of  the  cathode  porosity  values.  As  the 
anode  porosity  decreased  to  0.12,  the  predicted  polarization  curves 
demonstrated  a  severe  mass-transport  potential  loss  beyond 
a  current  density  of  100  mA  cm-2  (Fig.  4(b)).  The  predicted  Pmax 
was  42  mW  cm-2,  a  40%  reduction  from  the  value  obtained  for  an 
anode  porosity  of  0.41.  Again,  the  data  confirm  that  the  anode 
porosity  is  more  influential  than  the  cathode  porosity. 

4.2.  Comparison  between  the  proposed  and  conventional 
Kulikovsky  model 

Fig.  5  presents  the  DMFC  polarization  curves  using  the 
conventional  Kulikovsky  model  [38]  and  our  modified  model,  along 
with  experimental  data  for  comparison.  Fig.  5(a)  consists  of  the 
experimental  data  and  the  model  predictions  using  a  Nafion-117 
membrane  electrolyte  and  fed  with  1  M  methanol  as  the  anode 
feed  at  70  °C.  The  modified  model  fits  well  with  the  experimental 
data  reported  by  Ge  and  Liu  [39],  whereas  the  Kulikovsky  model 
suggests  much  lower  cell  voltage  and  power  density  values.  The 
Pmax  values  obtained  using  the  Kulikovsky  model  and  our  modified 
model  were  76  and  83.4  mW  cm-2,  respectively.  The  latter  was 
closer  to  the  experimental  value  of  86.7  mW  cm-2.  The  percentage 
errors  for  using  the  Kulikovsky  model  and  the  modified  model 
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Fig.  3.  Predicted  polarization  curves  using  the  modified  model  for  DMFCs  with 
a  porosity  of  (a)  0.12  and  (b)  0.41  for  the  cathode  diffusion  layer  and  various  porosities 
for  the  anode  diffusion  layer.  (The  cell  was  fed  with  1  M  methanol  and  pure  oxygen  at 
60  °C  and  was  equipped  with  a  Nafion-117  electrolyte.) 


Current  density  (mA  cm  2) 


Fig.  4.  Predicted  polarization  curves  using  the  modified  model  for  DMFC  with 
a  porosity  of  (a)  0.41  or  (b)  0.12  at  the  anode  diffusion  layer  and  various  porosities  at 
the  cathode  diffusion  layer.  (The  cell  was  fed  with  1  M  methanol  and  pure  oxygen  at 
60  °C  and  was  equipped  with  Nafion-117  electrolyte.) 


were  12.3  and  3.8%,  respectively.  Except  for  the  low  current  density 
regions,  the  modified  model  agrees  well  with  the  experimental 
data.  Our  modified  model  provides  precise  estimates  of  cell  voltage 
and  power  density  for  DMFCs. 

4.3.  Effect  of  temperature  on  cell  performance  using  the  modified 
model 

The  proposed  mathematical  model  for  DMFC  performance  was 
also  compared  with  the  experimental  data  at  80  °C  and  90  °C  using 
a  Nation  membrane  and  a  1  M  methanol  anode  feed.  The  experi¬ 
mental  results  were  obtained  from  Gurau  and  Smotkin  [40],  Scott 
et  al.  [41],  Wu  and  Zhang  [42],  and  Witham  et  al.  [43]  and  are 
plotted  in  Fig.  5(b)  and  (c).  The  voltage  and  power  density  cure 
pattern  of  the  modified  model  for  the  Nation  membrane  were  in 
agreement  with  the  experimental  values  of  Gurau  and  Smotkin  [40] 
and  Scott  et  al.  [41  ].  Again,  the  modified  model  outperformed  the 
conventional  model. 

Fig.  6(a)  and  (b)  plots  the  predicted  cell  voltage  and  power 
density  as  functions  of  current  density  using  the  modified  model 
for  DMFCs  equipped  with  Nation  at  40-90  °C  to  verify  the  accuracy 
of  the  modified  model.  As  the  cell  temperature  increased  from 
40  °C  to  90  °C,  Pmax  increased  from  48.8  to  112.14  mW  cm”2.  This 


increase  resulted  from  the  higher  Nation  conductivity  associated 
with  a  higher  temperature,  leading  to  a  slower  decline  in  ohmic  loss 
(Fig.  6(a)).  The  predicted  Pmax  data  using  the  modified  model  were 
in  excellent  agreement  with  the  experimental  results  for  the  entire 
temperature  range,  as  demonstrated  in  Fig.  6(c).  One  of  the  main 
contributions  of  this  approach  is  to  predict  DMFC  performance  for 
different  temperatures  as  long  as  the  diffusion  layer  porosity  data 
are  available.  Other  determinations  are  based  on  theoretical  mass- 
transfer  equations.  One  feasible  approach  is  to  estimate  the 
porosity  from  one  temperature  at  which  the  effective  diffusion 
coefficients  are  determined.  By  assuming  that  the  porosities  are  not 
significantly  affected  by  the  temperature  change,  the  porosity 
values  can  be  used  to  estimate  the  effective  diffusivity  at  other 
temperatures  using  the  procedure  outlined  in  Section  2.2. 

4.4.  Effect  of  methanol  concentration  on  cell  performance  using 
modified  model 

Fig.  7(a)  and  (b)  plots  the  polarization  curves  and  the  power 
densities  obtained  using  the  modified  model  for  DMFCs  with  Nation 
membrane  with  different  methanol  feed  concentrations  at  70  °C.  The 
measured  methanol  permeability  from  various  concentrations  was 
7.77  x  10“6, 8.15  x  1CT6, 8.52  x  1CT6,  and  1.25  x  1CT5  cm2  s*1  for  1. 2. 
3,  and  5  M,  respectively,  at  70  °C.  The  estimated  Pmax  values  for  1, 2, 3, 
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Fig.  5.  Conventional  and  modified  model  predictions  and  comparison  with  experi¬ 
mental  data  at  various  temperatures.  The  cell  was  fed  with  1  M  methanol  and  pure 
oxygen  and  was  equipped  with  a  Nafion-117  electrolyte,  (a)  70  °C,  experimental  data 
from  Ref.  [39];  (b)  80  °C,  experimental  data  from  Refs.  [40,41];  (c)  90  °C,  experimental 
data  from  Refs.  [42,43]. 


and  5  M  methanol  feeds  were  85, 99, 115,  and  110  mW  cnrT2  at  70  °C. 
Compared  with  the  experimental  data  reported  by  Ge  and  Liu 
[39]  under  comparable  conditions,  the  Pmax  values  were  87,  132, 
137,  and  128  mW  cnrT2,  respectively,  at  the  corresponding 
methanol  concentrations.  The  modified  model  predicted  enhanced 
cell  voltage  and  power  density  with  increasing  methanol  concen¬ 
tration,  although  the  concentration  effect  is  more  profound  in  the 
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Fig.  6.  Predicted  (a)  polarization  curves,  (b)  power  densities,  and  (c)  peak  power 
density  data  along  with  the  experimental  values  (cited  from  Refs.  [39-43]  and  our 
measurements)  of  DMFCs  at  30-90  °C  using  the  modified  model.  (The  cell  was  fed 
with  a  1  M  methanol  solution,  pure  oxygen,  and  Nafion-117). 


experimental  data.  On  average,  the  modified  model  predictions 
deviated  by  approximately  16%  from  the  data  presented  by  Ge  and 
Liu  [39].  Further  investigation  is  suggested  to  improve  the  model 
accuracy  by  incorporating  parameters  associated  with  the  high 
methanol  concentrations. 
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Fig.  7.  Predicted  (a)  polarization  curves  and  (b)  power  density  results  using  the 
modified  model  for  DMFCs  at  70  °C  and  fed  with  a  methanol  aqueous  solution  of 
various  concentrations  [39]. 


40 


H  30  o 
1 

H  20  ^ 

in 

c 

CD 
TJ 


-  10 


a> 

£ 

o 

CL 


Fig.  8.  Predicted  polarization  curves  (left  axis)  and  power  density  results  (right  axis) 
using  SPPEK  electrolyte  in  DMFC  (fed  with  2  M  or  3  M  methanol  solution,  oxygen  at 
cathode,  operated  at  60  °C). 
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Fig.  9.  Predicted  polarization  curves  (left  axis)  and  power  density  results  (right  axis) 
using  SPEEK  electrolyte  in  DMFC  (fed  with  2  M  methanol  solution,  oxygen  at  cathode, 
operated  at  80  °C).  Electrolyte  characters  and  cell  performance  data  are  from  Li  et  al. 
[44]. 


4.5.  Modeling  for  other  membrane  electrolytes 

The  DMFC  experimental  results  using  the  sulfonated  poly(- 
phthalazinone  ether  ketone)  (SPPEK)  were  used  to  verify  the 
versatility  of  the  modified  model.  This  electrolyte  was  prepared  and 
characterized  in  house.  Fig.  8  presents  the  DMFC  performance 
predictions  and  the  experimental  data  using  the  SPPEK  with  2  M 
and  3  M  methanol  feeds  at  60  °C.  The  estimated  Pmax  using  the 
SPPEK  was  26.2  mW  cm-2,  matching  the  experimental  result  of 
27  mW  cm-2  with  the  2  M  methanol  feed.  For  the  3  M  methanol 
feed,  the  Pmax  prediction  was  30.2  mW  cm-2,  in  agreement  with  the 
experimental  value  of  30.4  mW  cm-2.  The  predicted  peak  power 
density  results  for  both  the  2  M  and  3  M  methanol  feeds  are  almost 
identical  to  the  experimental  results.  This  agreement  confirms  the 
accuracy  and  versatility  of  the  modified  model. 

We  also  validated  our  proposed  model  using  data  presented  by 
Li  et  al.  [44].  These  researchers  prepared  and  characterized 
a  sulfonated  poly(ether  ether  ketone)  (SPEEK)  membrane,  followed 
by  determining  DMFC  performance  using  this  solid  electrolyte. 
Fig.  9  illustrates  the  prediction  and  experimental  results  for  the  cell 
voltage  and  power  density  for  the  DMFC  using  the  SPEEK  electro¬ 
lyte.  The  peak  power  density  prediction  was  69  mW  cnrT2,  slightly 
higher  than  the  experimental  value  (64  mW  citT2).  However,  the 


error  was  less  than  7.3%.  This  example  demonstrates  that  given  the 
characteristics  of  the  electrolyte  (thickness,  conductivity,  and 
methanol  permeability),  one  can  apply  our  modified  model  to 
estimate  the  DMFC  performance.  This  modeling  process  provides 
cost-saving  advantages  by  eliminating  unnecessary  fuel  cell  prep¬ 
aration  and  testing  experiments,  particularly  by  reducing  the  usage 
of  expensive  catalysts  and  electrolytes.  Preliminary  data  on  fuel  cell 
performance  can  be  simulated  with  sufficient  precision  using  this 
model. 

5.  Conclusion 

The  effective  diffusion  coefficients  of  reactants  in  the  porous 
diffusion  layers  were  estimated  using  the  Wilke-Chang  equation, 
the  Lennard-Jones  potential  model,  and  the  UNIQUAC  equation  for 
vapor-liquid  equilibrium  data  on  methanol/water  solutions.  The 
data  were  adopted  into  a  conventional  Kulikovsky  model  to  esti¬ 
mate  the  polarization  curve  and  power  density  of  a  DMFC.  There¬ 
fore,  the  prediction  has  a  theoretical  basis  without  the  requirement 
of  any  fitting  parameters.  The  model  was  validated  using  experi¬ 
mental  results  obtained  using  a  1-5  M  methanol  feed,  operating 
temperatures  at  30-80  °C,  with  Nation  and  other  solid  electrolytes. 
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This  model  was  able  to  predict  the  effects  of  diffusion  layer  porosity, 
cell  temperature,  and  methanol  concentration  on  cell  performance 
using  various  electrolyte  membranes.  The  data  in  Figs.  2  and  5-9 
illustrate  the  precision  of  the  predictions  (shown  in  lines)  against 
experimental  results  from  various  sources  (shown  in  symbols).  The 
simulation  model  provides  a  valuable  development  tool  with  cost¬ 
and  time-saving  advantages  for  researchers. 
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Nomenclature 

a:  activity  of  reactant  or  product  (— ) 
ba:  Tafel  slope  at  anode  (V) 
bc:  Tafel  slope  at  cathode  (V) 

Cm:  methanol  concentration  (mol  cm-3) 

Co:  oxygen  concentration  (mol  cm-3) 

D°:  solute  diffusion  coefficient  at  infinite  dilution  (cm2  s-1) 

Dba:  methanol  diffusion  coefficient  in  anode  diffusion  layer  (cm2  s-1) 

Dt,c:  oxygen  diffusion  coefficient  in  cathode  diffusion  layer  (cm2  s-1) 

Djj:  diffusion  coefficient  of  component  I  in  mixture  with  J  (cm2  s-1 ) 

D6^:  effective  diffusion  coefficient  of  reactant  in  porous  diffusion  layer  (cm2  s-1) 

E:  thermodynamic  cell  potential  at  any  condition  (V) 

E°:  thermodynamic  cell  potential  at  standard  state  (V) 

F:  Faraday’s  constant  (96,485  C  mol-1) 

AG°:  Gibbs  free  energy  change  in  overall  reaction  for  a  DMFC  at  standard  state 
(J  mol-1) 

i:  current  density  (A  cm-2) 
iea:  exchange  current  density  for  anode  (A  cm-2) 
iec:  exchange  current  density  for  cathode  (A  cm-2) 
i;a:  limiting  current  density  at  the  anode  (A  cm-2) 
ijc;  limiting  current  density  at  the  cathode  (A  cm-2) 

Lba '■  thickness  of  backing  layer  at  the  anode  (cm) 

LjjC:  thickness  of  backing  layer  at  the  cathode  (cm) 

L:  thickness  of  electrolyte  (cm) 

M:  molar  mass  (g  mol-1) 
n;  number  of  electrons  in  half-reaction 
P:  power  density  (mW  cm-2) 
p;  gas  pressure  (bar) 

R:  gas  constant  (8.314  J  mol-1  K-1) 

Rc:  dimensionless  parameter  accounting  for  potential  drop  due  to  methanol  cross¬ 
over  (— ) 

T:  temperature  (K) 

T*:  reduced  temperature  (— ) 

V:  cell  voltage  (V) 
v:  molar  volume  (cm3) 
x:  mole  fraction  (— ) 

Greek 

a:  thermodynamic  correction  term  (— ) 

(3:  methanol  permeability  (cm2  s-1) 
e:  characteristic  Lennard-Jones  energy  (J) 
solvent  viscosity  (cP) 
r\a:  overpotential  at  the  anode  (V) 
rjc:  overpotential  at  the  cathode  (V) 
t)q:  ohmic  potential  drop  (V) 

6:  porosity  of  diffusion  layer  (— ) 
k:  Boltzmann’s  constant  (J  K-1) 
p:  dimensionless  parameter  (— ) 

a:  proton  conductivity  of  solid  electrolyte  membrane  (S  m-1) 
goa:  characteristic  Lennard-Jones  length  (A) 

0:  solvent  association  constant  (— ) 

Qd:  diffusion  collision  integral  (— ) 

Subscript 
a:  anode 
A:  air 
c:  cathode 
M:  methanol 
O;  oxygen 
ref:  reference 
W:  water 


